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Abstract 

We present findings for DG Tau B and W CrA, two of the objects observed in our Spitzer IRS 

project to search for molecular absorption in edge-on disks, along with near-IR spectroscopy of the 
CO fundamental transitions and mid-IR imaging. While the only gas absorption seen in the Spitzer 
IRS spectrum toward DG Tau B is CO2, we use gas abundances and gas/ice ratios to argue that 
we are probing regions of the disk that have low organic molecule abundances. This implies the 
rarity of detecting molecular absorption toward even edge-on disks with Spitzer IRS is a result of 
high dependence on the line of sight. We also argue the disk around DG Tau B shows high amounts 
of grain growth and settling. For VV CrA, we use the silicate absorption feature to estimate a dust 
extinction, and model the disk with a spectral energy distribution fitting tool to give evidence in 
support of the disk geometry presented by Smith et al. (2009) where the Primary disk is the main 
source of extinction toward the infrared companion. 

Subject headings: circumstellar matter - ISM: abundances - planetary systems: protoplanetary disks 
- stars: individual (DG Tau B, VV CrA) - stars: pre-main-sequence 



1. INTRODUCTION 

The study of disks of dust and gas surrounding young 
stars provides insight into the environment in which plan- 
ets are likely to be forming. Molecular abundances in the 
disk can be used to understand the chemical evolution of 
the gas, disk processes such as radial and vertical mixing, 
and the effects of stellar irradiation on the disk chemistry 
(Bergin et al. 2007). While models predict disk chem- 
istry in the earth-like planet forming regions (i.e. Mark- 
wick et al. 2002; Agundez et al. 2008; Glassgold et al. 
2009; Woods & Willacy 2009), supporting observations 
of the chemical tracers have been scarce to date. 

The Spitzer Infrared Spectrograph (IRS) made obser- 
vations that probed the chemistry of the inner disk, de- 
tecting molecular absorption (Lahuis et al. 2006) and 
emission (Carr & Najita 2008; Salyk et al. 2008). The 
resolution of IRS was well-suited for detecting molecules 
with large equivalent widths from ro-vibrational Q- 
brariches. Lahuis et al. (2006) detected warm C2H2, 
HCN. and CO2 in absorption toward IRS 46, likely orig- 
inating in the inner regions of a circumstellar disk seen 
near edge-on. Unfortunately, systems displaying molec- 
ular absorption are rare with, for example, IRS 46 being 
the only one of roughly 100 young stellar objects (YSO) 
observed in the c2d Spitzer Legacy program (Evans et 
al. 2003). 

Molecular C2H2 and HCN absorption have also been 
detected toward GV Tau N (Doppmann et al. 2008; cf. 



Gibb et al. 2007), a classical infrared companion (IRC; 
Koresko et al. 1997). IRCs are binary companions to 
young stars that are faint in the visible but dominate 
the system flux in the infrared. Correia et al. (2007) 
found with MIDI/VLTI data that three IRCs, including 
GV Tau, are likely T Tauri stars (TTS) with the disk 
seen almost edge-on. Doppmann et al. (2008) found the 
warm temperature and the relative velocity shifts of the 
molecular absorption further supported the edge-on disk 
model. Gibb et al. (2007) compared the molecular ab- 
sorption abundances to predictions made by the chemical 
disk model of Markwick et al. (2002) and found them to 
be consistent with the inner region of a protoplanetary 
disk. 

Because IRS 46 and GV Tau N show the rare features 

that are valuable for probing the inner disk chemistry, 
it is advantageous to study more disks with similar at- 
tributes. We were granted time with Spitzer IRS to ob- 
tain high signal-to-noise (S/N) spectra of seven YSOs. 
Our target selection was limited to systems thought to 
be similar to IRS 46 and GV Tau N with a large column 
of gas between the observer and the emitting source nec- 
essary for absorption studies: three circTinistellar disks 
seen near edge-on and four classical IRCs. DG Tau B and 
W CrA are two of the objects observed in this project. 

Mundt & Fried (1983) first designated DG Tau B as 
the source of a Herbig-Haro jet (HH 159) that flows at 
15° to the plane of the sky. Near-infrared HST images 
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show DG Tau B to be a bipolar reflection nebula em- 
bedded in an envelope (Stapelfeldt et al. 1997). Stark 
et al. (2006) modeled near-infrared disk images, flnding 
it likely has a disk inclination of ^73°. Previous Spitzer 
IRS measurements have shown this system displays CO2 
ice (Watson et al. 2004) and gas absorption (Pontopp- 
idan et al. 2008), so we obtained higher S/N to search 
for other species. 

The northeast component of the VV CrA system is a 
classical IRC, and Correia et al. (2007) used MIDI/VLTI 
observations to show the IRC is likely a TTS seen 
near edge-on. Smith et al. (2009; hereafter S09) used 
CRIRES/'VLT to detect CO gas absorption with mul- 
tiple temperature components toward the Primary and 
IRC. They proposed two disk morphology scenarios to 
explain the CO absorption seen toward VV CrA IRC: 
(A) the gas absorption is in a disk around the Primary 
in the line of sight to the IRC or (B) all the CO absorp- 
tion is in a disk around the IRC. While S09 considers 
Case A to be more likely, with the observed CO residing 
in the warm surface layer and cooler midplane of the Pri- 
mary disk, neither case could be dismissed so we consider 
both geometries. 

In addition to the IRS spectra, we observed both tar- 
gets with ground-based telescopes. We used M-band 
spectroscopy to observe the CO fundamental (Av=l) 
transitions toward both objects. Observing these lines 
probes a wide range of radii in circumstellar disks and 
has proven useful for exploring the disk gas in the inner 
disk regions (Najita et al. 2003; Brittain et al. 2007 
and references therein). Finally, because molecular ab- 
sorption is limited to studies in the line of sight to the 
continuum emitting source and the Spitzer IRS spectrum 
combines both VV CrA components, we obtained mid- 
IR images of W CrA in order to map the continuum 
at spatial scales unavailable with Spitzer and to separate 
the continuum flux for the two components. 

We present here results for DG Tau B and VV CrA 
from our study of sources in the Spitzer survey. We in- 
clude high-resolution absorption spectroscopy and, for 
W CrA, mid-IR imaging. In §2 we introduce the obser- 
vations taken; in §3 we describe the data reduction and 
discuss the molecular detections; and we discuss model- 
ing and the conclusions in §4. 

2. OBSERVATIONS 

2.1. Spitzer IRS 

We observed DG Tau B and VV CrA as part of our Cy- 
cle 5 General Observers project. Program ID 50152. We 
used the Spitzer Infrared Spectrograph (IRS; Houck et 
al. 2004) short-low and short-high modules (SL and SH) 
in stare mode to search for molecular absorption species 
with large equivalent widths such as ro-vibrational Q- 
branches. The SL spectra obtained have a coverage 5.2- 
14.5 \xm with R=A/AA=60-127, and the SH spectra have 
a coverage of 9.9-19.6 ^m with R=600. We adopted the 
observing technique described in Carr & Najita (2008). 
In order to verify the repeatability of the observations, 
we lowered the exposure time for each data frame and 
increased the number of integrations. Each setting (SLl, 
SL2, SH) thus had twelve on-source data frames at each 
nod position, each with a 6 second integration. For hot 
pixel subtraction and to remove background emission. 



each observation was supplemented by off-source inte- 
grations for half the time as on-source (Carr & Najita 
2008). DG Tau B (AOR: 25680128) was observed on 
2008 November 5 at 4''27™2.66'* +26°5'30.5" (J2000). 
VV CrA (AOR: 25680896) was observed on 2008 Novem- 
ber 18 with the pointing on the IRC at 19'*3"6.86-' - 
37°12'48.2" (J2000). The coordinates for VV CrA were 
chosen to center the slit on the IRC, and were calculated 
using a 2.1" offset, 44° east of north (CheUi et al. 1995) 
from the known primary coordinates found on SIMBAD. 
Both observations used the nominal pointing accuracy of 
0.1". 

2.2. Infrared Spectroscopy 

On 2008 November 16, we used NIRSPEC (McLean 
et al. 1998) on Keck II in echelle mode with R=25,000 
to detect the CO fundamental transitions toward DG 
Tau B. Our instrument setup included the M-wide block- 
ing filter, along with echelle and cross disperser settings 
to image orders 14-16 so we could view the ^^CO v=0- 
1 R(0)-R(9) and P(25)-P(30) transitions. We used the 
standard ABBA nod sequence with thirty 1-second ex- 
posures co-added at each position. Standard star obser- 
vations were taken before and after DG Tau B for telluric 
division. 

For a similar observation of both VV CrA components, 

we used Phoenix (Hinkle et al. 1998) on Gemini South 
on 2009 April 6. We used the standard four pixel slit 
(0.35"), giving R=50,000. The M2150 filter and instru- 
ment settings were used to detect the ^^CO v=l-0 R(3)- 
R(5) transitions near 4.66 ^m. Each ABBA nod position 
had a 120 second integration (four 30 second exposures 
co-added internally), and we used a telescope position 
angle of 44° to place both binary components within the 
slit. Standard star observations were taken after the sci- 
ence frames for telluric division. 

2.3. Imaging 

We imaged the binary components using T-ReCS 
(Telesco et al. 1998) on Gemini South using 6 filters rang- 
ing from 7.7-25.5 [im (see Table 1). We cycled through 
each of the filters while staying on the science target, 
and took calibration images of standard stars before and 
after the science target for point spread function (PSF) 
corrections and to test for any changes in PSF or cali- 
bration during the target observation. Observations were 
conducted using the standard chop-nod mode. When we 
nodded, the target was placed at the center of the field 
of view for half the total integration time, and the other 
half placed the target into a corner of the array for frames 
to be used for background subtraction. 

3. DATA REDUCTION AND DETECTIONS 

3.1. Spitzer IRS 

We used the reduction technique described in Carr &: 
Najita (2008). We used the integrations of the nearby 
blank sky to identify hot pixels, which were removed 
from the science frames by interpolation using custom 
routines in IDL. We then used IRAF to extract the raw 
spectra from each individual raw data frame. A similar 
procedure was used to extract spectra from the regularly 
observed HR 6688, which we used for standard star di- 
vision, to remove fringing, and for fiux calibration. Be- 
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cause the fringing in each observation is dependent on 
the target position on the sht, each data frame was di- 
vided by the standard star observation that showed the 
least amount of fringing in the shorter wavelengths after 
division. Any residual fringing was removed using IRS- 
FRINGE in IDL. Calibration of the SL spectra is highly 
dependent on the slit position, so we used the optimal 
point source extraction in the Spitzer IRS Custom Ex- 
tractor (SPICE; version 2.2) to extract the individual 
spectra before co-adding. 

The Spitzer IRS spectrum of DG Tau B, with an aver- 
age S/N=100 per pixel, can be seen in Fig. 1. A 10% con- 
tinuum flux difference between the SH and SL modules 
was corrected by scaling the SH by a constant, without 
any adjustment to the continuum shape, to match the 
SL at wavelengths where the orders overlap. The spec- 
trum shows the amorphous silicates, CO2 ice absorption, 
water ice, and an absorption feature at 6.8 jim, which is 
possibly due to methanol (Watson et al. 2004, Zasowski 
et al. 2009), as well as the CO2 gas absorption at 15.2 ^m 
as previously recognized by Pontoppidan et al. (2008). 
We further detect [Nell] emission at 12.81 [Lm with line 
intensity 3.0 ± 1.5 x 10"^^ ergs cm~^ s~^. 

The Spitzer IRS spectrum of VV CrA shown in Fig. 
1, with S/N=160 per pixel, is of the combined binary 
components as they were unresolved within the sHt. It 
can be seen there is a continuum mismatch between SH 
and SL modules, likely due to being centered on the IRC 
and from the modules having different slit position an- 
gles and widths. The SH and SL modules have posi- 
tion angles 221.5° and 136.7° (for a Spitzer roll angle 
of 0°), respectively, so during the SH observation the 
binary components were separated 0.31" across the 4.7" 
slit while during the SL observation the components were 
separated 1.20" across the 3.6" slit. The SLl is at an 
intermediate wavelength and flux, so we scaled the SL2 
and SH modules by a constant to match the intermediate 
SLl continuum. We find the system displays amorphous 
silicate absorption and CO2 ice absorption at 15.2 ^m. 

3.2. CO Fundamental 

To reduce our Phoenix and NIRSPEC spectra, we used 
standard spectral reduction routines in IRAF. After flat- 
fielding each frame, we subtracted nod pairs to remove 
the sky emission. The target spectra were extracted us- 
ing apall in IRAF with background subtraction using 
both sides of the profile (from both sides of the total 
binary profile for VV CrA). In order to increase S/N 
and for better telluric division, we first extracted the 
spectrum from each individual science and standard star 
frame (standards listed in Table 1). We then divided 
each science spectrum by each standard star spectrum 
using their corresponding airmass, then co-added the 
spectra. We found that applying telluric division before 
co-adding the spectra resulted in better telluric correc- 
tion. The wavelength solution was made by using the 
IRAF task identify on sky emission lines in coadded sci- 
ence frames. The rest frame wavelengths of the emission 
lines were taken from the HITRAN database (Rothman 
et al. 1998). 

Our NIRSPEC spectra of DG Tau B, seen in Fig. 2, 
have S/N=21 and 28 per pixel for orders 15 and 16, re- 
spectively, and we found both ^^CO and ^^CO in ab- 



sorption (see Table 2). The average line profiles for the 
low- and high-J ^^CO transitions have slightly different 
velocity shifts and widths (see Table 3 and Fig. 3) which 
indicates we are likely probing multiple velocity compo- 
nents. 

We were able to clearly resolve the VV CrA compo- 
nents with Phoenix, and the spectra, shown in Fig. 4, 
had an average S/N of 220 per pixel for the IRC and 
280 for the Primary. The W CrA IRC displays ^^CO 
and optically thick -'^^CO in absorption (see Table 4), as 
well as a shallow, broad ^^CO absorption feature with 
VLSfl = -24 km s"^ and FWHM=32.9 km s~^ (see line 
profile in Fig. 5). The VV CrA Primary shows cold ^^CO 
within broader ^^CO emission lines. The cold, unresolved 
^^CO absorption hues have FWHM=5.35 km s^^ while 
using the 4-pixel slitwidth with resolution R=6 km s~^, 
suggesting our resolution was limited by the PSF of the 
target rather than the entire slit width. 

The Primary and IRC both also show a complicated 
emission structure. This structure can be understood 
as a blend of broad CO emission lines from numerous 
ro-vibrational levels. Both of the W CrA components 
show CO overtone emission (Prato et al. 2003) with 
emission up to at least the v=4-2 vibrational band. The 
observed CO fundamental bands are consistent with the 
expectation of a hot, optically thick emission spectrum. 
In this respect, the spectra are similar to the extreme T 
Tauri star RW Aur (Najita et al. 2003). 

3.3. Imaging 

The T-ReCS images show a resolved VV CrA binary 
(Fig. 6). Images were made by subtracting nod pairs 
and co-adding the resulting images. High winds caused 
poor image quality, so we created custom IDL routines to 
center each data frame on the peak of the primary com- 
panion before co-adding. This gave better final image 
quality, but we did not detect any extended emission. 
We measured photometric fluxes using a 2.1" aperture 
centered on each component and on the standard star 
observations, which were then used for flux calibration 
(see results in Table 5). While the primary does show a 
hint of silicate absorption, a majority of this feature seen 
in the Spitzer IRS spectrum originates in the IRC. We 
flnd that the broadband fluxes measured with T-ReCS fit 
the Spitzer IRS spectrum if we assume we are collecting 
only 70% of the flux from the Primary (see Fig. 1), with 
the flux loss likely due to slit losses from being centered 
on the IRC. 

4. DISCUSSION 

4.1. Molecular Detections 

4.1.1. DG Tau B 

All of our models of the CO and CO2 gas were done 
using a single temperature slab model. We fit the ^^CO 

R(16-19) transitions with the assumption the transitions 
are optically thin and found a temperature of 288±7 K 
with column density N(^'^CO)=3.5±G.6 x 10^'' cm^^ (see 
population diagram in Fig. 7). If we fit the high-J ^^CO 
P(25-27, 30) lines to an optically thin model as well, we 
find the gas would be at a high temperature of 930 ± 
60 K with N(^^CO)=7.86 ± 0.2 x 10^^ cm'^. However, 
this would indicate a ^^CO/^^CO ratio of only 2.2, much 
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lower than the mean ratio found in the local ISM (Wilson 
1999). This indicates all of the observed ^^CO lines arc 
likely optically thick. To constrain the temperature of 
the absorbing gas, we modeled ^^CO gas with a column 
density 2.42 x 10^^ cm'^^ estimated using the N(i-'^CO) 
and mean isotopic ratio of 69±6 from Wilson (1999), and 
varied the turbulence and temperature. We could not 
model all the ^^CO transitions to a single temperature, 
with the low-J transitions always showing higher optical 
depths than permitted by fits to the high-J transitions. 
This is likely due to multiple velocity components being 
viewed in the low-J transitions, which would account for 
the difference in velocity shift and width from the high- 
J transitions. We thus fit only the high-J P (25-27, 30) 
transitions, and the goodness-of-fit, shown in Fig. 7, 
indicates the gas is likely 370 ± 20 K with a Doppler 
b- value vt„rt = 2.1 ± 0.2 km s~\ If the ^^CO is self- 
shielding and the ^^CO/'^'^CO is raised to 100, similar to 
the ratios found in S09, these values would not change 
much as the temperature would still be ~ 350 K and 
vt„rb ~ 2 km s~^ 

We find the model for the CO2 gas has degenerate fits 
by changing the microturbulence and temperature, so we 
use the Doppler b- value of 2.1 km s^^ found with the CO 
gas, and find T=300 K and N(C02) = 2.9 x 10^'^ cm.-^. 
We further find upper limits for C2H2 and HCN, two 
other organic molecules detected toward IRS 46 (Lahuis 
et al. 2006), to compare molecular abundance ratios. 
To get upper limits on column densities, we note that 
temperatures of both the CO and CO2 gases toward DC 
Tau B are similar to that found toward IRS 46 (400 and 
300 K, respectively), so wc start with the assTimption 
that if C2II2 and HCN are present, they would have simi- 
lar temperatures to IRS 46 (700 and 400 K, respectively). 
To derive the upper limits, wc individually modeled the 
absorption spectrum for the molecules at those temper- 
atures. We then used the model-generated profiles to 
measure feature strengths at various wavelengths in the 
DC Tau B spectrum, giving us a distribution of feature 
strengths we would expect to measure given only the 
spectral noise. Wc then set the upper limit as the col- 
umn density that would provide a 3a integrated intensity. 
We find the upper limit on the C2H2/CO2 ratio (< 0.19) 
is not very different from the ratio measured for IRS 46 
(0.3), suggesting the non-detection of C2II2 in our data 
may be the result of the lower absorbing column in this 
source. The upper hmit on the HCN/CO2 ratio (< 0.22) 
is a factor of two smaller than was measured for IRS 46 
(0.5). Modeling the gas at higher temperatures with sim- 
ilar column densities decreases the integrated intensity 
of the absorption features, so higher gas temperatures 
would result in higher upper limits on the cohimn den- 
sities. A temperature increase of 50 K would raise the 
column density upper limit by ~15% for both molecules. 
If we model the molecules at 300 K, like the CO2 gas, 
the upper limits of the C2H2/CO2 and HCN/CO2 ratios 
would be 0.06 and 0.18, respectively. 

The [Nell] emission, as seen toward DC Tau B, is a 
feature commonly found by Spitzer IRS toward TTSs 
(Lahuis et al. 2007). Giidel et al. (2010) showed that 
among TTS the [Nell] emission strength increases with 
the jet mass loss rate, the latter as inferred from the high 
velocity component of the [01] emission. Mundt & Fried 



(1983) found that HH 159, a jet associated with DC Tau 
B, has a mass loss rate of 10~^ Mq yr^^. If we assume 
DC Tau B is 145 pc away, our measured [Nell] luminosity 
of 7.5 X 10^^ ergs s~^ is similar to that of TTS with com- 
parable mass loss rates (cf. Fig. 5, Giidel et al. 2010). 
As we do not see C2II2 and HCN in the disk, a simple ex- 
planation is they have low abundances in the disk. This 
could possibly be due to photodissociation. The [Nell] 
emission may be an indication the gas is irradiated by 
X-rays (Glassgold et al. 2007), and/or UV irradiation 
(HoUenbach & Gorti 2009), which in turn could reduce 
the C2H2 and HCN in the disk. However, many T Tauri 
stars with [Nell] emission similar to or stronger than that 
of DG Tau B show emission in C2H2 and HCN (Carr & 
Najita 2011, in prep). 

Alternatively, as C2H2 and HCN are commonly found 
in circumstellar disks (Pontoppidan et al. 2010), these 
species could be in the disk but simply not in our line 
of sight. This may be implied by the most striking dif- 
ference between DG Tau B and IRS 46: while the col- 
umn density of CO toward DG Tau B is only 10 times 
that toward IRS 46, the gas ratios to CO in DG Tau 
B are l%-2% of that found toward IRS 46 (see Table 
6). Also, the limits on the ratios N(C2H2)/N(CO) and 
N(HCN)/N(CO) are much smaller for DG Tau B than 
the measured ratios toward GV Tau N (Gibb et al. 2007). 
While the large CO column toward DG Tau B implies we 
are probing a large amount of disk material, we are not 
detecting the large columns of other gaseous species. 

The disk chemistry model in Agundez et al. (2008) 
may provide further insight into the disk conditions 
probed by our observations. They focus their models 
on the photon dominated region of the disk, and use 
photochemistry to estimate the steady state abundances 
of simple organic molecules. They show the C2H2 and 
HCN have low abundances for disk radii >1.5 AU, and 
the CO/CO2 is large for radii <1 AU. While the tem- 
peratures we measured for CO and CO2 are cooler than 
expected for radii <1 AU, the vertical structure of these 
molecules is not available. One possible scenario is the 
line of sight passes below the main temperature inversion 
layer. 

Doppmann et al. (2008) note GV Tau N has a larger 
TSio/Tco^ice than other Class I YSOs in Taurus (typi- 
cally < 4; Furlan et al. 2008), and suggest this could 
happen if a majority of the probed silicate is residing 
in a warm region of the disk atmosphere. We calculate 
this ratio for DG Tau B and IRS 46 for comparison. To 
do this, we used a polynomial fit across the features to 
determine the initial continua. For the silicate feature, 
we fit to the continuum in the 13-14.8 ^m region, and 
to the continuum between the water and "methanol" ice 
at wavelengths < 7.8 [ivn. For the CO2 feature, we as- 
sume the CO2 ice is located at larger disk radii than the 
silicate, as found by Watson et al. (2004), so the ice is 
absorbing from a continuum already including the sili- 
cate absorption features. We thus use a continuum for 
DG Tau B that is modified from that used in Watson et 
al. (2004), found by fitting the spectrum at wavelengths 
13-14.8 (im and >15.6 \mi. The continuum fits are shown 
in Fig. 1. We then used the relation l^bs = IeTOexp(— r), 
where Io6s is the observed flux and lem is the emitted con- 
tinuum, to flnd Tsio = 1.52±0.06 and rcos = 0.24±0.01. 
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Lahuis et al. (2006) found IRS 46 to have tsw = 0.86, 
and Pontoppidan et al. (2008) shows tco-j =0.1. While 
both IRS 46 and DG Tau B have higher ratios than found 
for Class I objects in Taurus, both GV Tau N and IRS 
46 have higher ratios (~ 9 and 8.6, respectively) than 
DG Tau B (6.6±0.4). Additionally, using the CO2 ice 
column density measurements for IRS 46 and DG Tau B 
found in Pontoppidan et al. (2008), we find IRS 46 has 
a higher CO2 gas/ice ratio, about 8 times that of DG 
Tau B. This suggests the line of sight toward IRS 46 is 
predominantly probing warmer regions than toward DG 
Tau B, resulting in the detection of C2H2 and HCN. 

Alternatively (or in addition), the ratios of CO2 gas, 
ice and silicate could simply be a factor of disk incli- 
nation. If the lines of sight for IRS 46 and GV Tau N 
pass higher in the disk atmosphere than for DG Tau B, 
raising the silicate to ice ratio, it could also be probing 
a moderately different region of the disk where the or- 
ganic molecules reside. If we are only viewing the cold to 
moderate temperature (< 400 K) regions deeper in the 
disk in DG Tau B, the C2H2 and HCN production would 
be lower in these regions, thus lowering the steady state 
abundances (Agundez et al. 2008) and resulting in their 
non-detection. If this is the case, the detection of pre- 
biotic molecules in absorption may be highly sensitive to 
the line of sight, even for edge-on disks. This could pos- 
sibly further explain the rarity of molecular absorption 
among Spitzer IRS spectra of young stars. 

4.1.2. VV CrA 

Our Phoenix spectra of the W CrA system shows 
^^CO transitions in absorption toward both components. 
However, the optically thick lines for the IRC and few 
lines for the primary make it difficult to extract accu- 
rate temperature and column density measurements. We 
therefore use the CO abundances reported in S09 for the 
IRC. To get the temperature and column density of the 
cold i^CO toward the primary, we obtained the CRIRES 
spectra used in S09 from the ESO archive. We reduced 
the data using a similar procedure to that used for the 
Phoenix and NIRSPEC, although the wavelength cali- 
bration and telluric division were performed using a tel- 
luric model fitter (Seifahrt et al. 2010). We find the 
cold gas toward the primary has a column density of 
2.6 ± 0.4 X 10^^ cm~^ and temperature 29±4 K (popula- 
tion diagram is shown in Fig. 8). While doppler velocity 
shifts of the ^^CO v=(2-0) absorption lines reported for 
the IRC in S09 are ~ 3-5 km s"S the CRIRES spectrum 
shows a smaller shift for the v=(l-0) lines (^ 1 km s"'^), 
similar to the shift we found with the Phoenix spectrum 
(0.78 ± 0.5 km s^^). One change between the Phoenix 
and CRIRES spectra is the CRIRES spectrum does not 
show the broad, shallow lines, but rather displays other 
^^CO lines in the -40 to -80 km s^^ range (see Fig. 9). 
Thus, these broad lines are in a turbulent, varying region, 
possibly located in an outflow. 

While the near-IR spectrum of VV CrA is quite com- 
plicated, the Spitzer IRS spectrum of the W CrA binary 
shows only amorphous silicate and CO2 ice absorption. 
To determine the optical depths of the silicate and CO2 
ice, we used the broadband fluxes found with the T-ReCS 
images to estimate continua for the two components. The 
T-ReCS images indicate almost all of the silicate absorp- 



tion is toward the IRC, so this feature is not associated 
with the ISM. To get the optical depth of the silicate 
and CO2 ice, we flrst subtracted the Primary continuum 
from the Spitzer IRS spectrum, then we used the same 
technique used to get the optical depths for DG Tau B. 
The peak optical depth for the 9.7 y.m feature seen to- 
ward the IRC is tsw = 1-35 ± 0.04. Using the relation 
Av/rsto = 16.6 ± 2.1 from Rieke & Lebofsky (1985) 
for grains in the ISM, we estimate the extinction to be 
Ay = 22.4 ± 3 mag, similar to the estimate by Koresko 
et al. (1997) of 19 mag. 

For the CO2 ice feature, if all of the absorption for 
this feature is associated with the IRC, the optical depth 
would be TCO2 = 0.055 ±0.004, which would give a ratio 
TSio/Tco-2 — 24.8 ± 2. At flrst glance, this high ratio 
could imply an edge-on disk geometry when comparing 
to the Class I YSOs in Taurus (Furlan et al. 2008). How- 
ever, the Tsio/'Tco2 ratio found for the IRC is typical of 
values found by Gibb et al. (2004) for embedded, mas- 
sive YSOs, and found by Alexander et al. (2003) for low- 
and intermediate-mass YSOs. Further, the selection in 
Alexandar et al. (2003) included YSOs in the CrA cloud, 
so an edge-on geometry may not be necessary for the high 
Tsio/Tco2 ratio seen toward the IRC. 

4.2. SED Model Fits 

We want to characterize the disk structures of DG Tau 
B and VV CrA in order to have a context for the gas ab- 
sorption, so we used the online SED fltting tool described 
in Robitaille et al. (2006) to get estimates on the disk 
structure. This tool flnds a best flt to an input spectral 
energy distribution (SED) using a database of 200,000 
pre-computed radiation transfer models of circumstellar 
disks, with 14 physical parameters varied at 10 differ- 
ent disk viewing angles. This modeler also estimates the 
interstellar and circumstellar extinction separately. The 
user is able to input spectral points, along with the aper- 
tures used in the measurements, at wavelengths ranging 
from the UV to radio. The user can also use upper and 
lower limits with a "confldence" percentage that rates 
how much to penalize models that exceed the limit. The 
output disk models can be limited to those within a dis- 
tance and interstellar extinction range. The models are 
returned in order of increasing value, and while a 
cut-off is arbitrary and impossible to make formal when 
fltting well behaved, symmetric models to less behaved 
YSOs (Robitaille et al. 2007), we flnd for our disks the 
SED modeler provides good flts up to xLst ~ ^ ION. 

While Stark et al. (2006) created synthetic images to 
model the envelope and disk structure seen in HST near- 
IR images, giving them the advantage of directly mod- 
eling images that show disk structure, we can compare 
and contrast to the model estimates from the SED fltter. 
This can provide insight into the model results for VV 
CrA. If we compare the model extinction for W CrA 
to our extinction estimates using the optical depth Tsio-, 
we may be able to constrain the location of the dust, 
whether in the disk or in the interstellar medium, and 
test the two disk scenarios presented in S09. 

4.2.1. DG Tau B 

We constrain the models for DG Tau B to be between 
120 and 170 pc away (GottUeb & Upson 1969; Elias 1978; 
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Kenyon et al. 1994) with interstellar extinction up to 
Ay = 15 mag (Whittet et al. 2001). Our Spitzer IRS 
spectrum was included by inputting the flux at every 1 
^m with a 10% error. The SED model fits can be seen in 
Fig. 10, and the parameter results for the best fit model 
are listed in Table 7. 

The models predict a disk that is not seen edge-on 
(18°), which is not surprising as Robitaillc et al. (2007) 
found this SED fitting tool does not determine the disk 
inclination well for most models. Robitaille et al. (2007) 
used the SED fitting tool to investigate DG Tau B and 
found an edge-on orientation (inclination=87°) for those 
models with no envelope accretion. However, when enve- 
lope accretion - as has been found for DG Tau B by Stark 
et al. (2006) - was added, all inclinations were allowed. 
We thus exclude disk inclinations from the reported disk 
parameters. 

Stark et al. (2006) modeled DG Tau B with a cen- 
tral star with = 0.5 Mq: T* = 3800 K, and min- 
imum disk radius 0.03 AU. They calculated an enve- 
lope infall rate of 5.0 x 10~^ Mq yr~^, which would be 
2.2 X 10~^ Mq yr~^ after correcting for our model stel- 
lar mass of 0.1 Mq, similar to the SED model prediction 
of 1.3 X 10~^ Mq yr~-^. However, the cavity opening an- 
gle and disk mass for our models are consistently smaller 
than those found by Stark et al. (2006). 

If we again use the ratio Ay/rsio = 16-6 ± 2.1 from 
Rieke & Lebofsky (1985), we can use our measured op- 
tical depth Tsio = 1.52 ± 0.06 to estimate the extinc- 
tion toward DG Tau B to be Ay = 25.2 ± 3.3 mag. 
Watson et al. (2004) also used this absorption feature 
to estimate an extinction of 30 mag, which they as- 
sociate with the circumstellar disk. The best fit SED 
model shows a somewhat higher circumstellar extinc- 
tion of Ay = 36.1 mag, but the circumstellar extinc- 
tion ranges from 26.6 to 42 mag for the disks within 
the cut-ofi^, comparable to our extinction estimate of 
25.2 ± 3.3 mag. 

While the model results for cavity opening angle and 
disk mass are small, the other parameters appear to 
be comparable to previous results and estimates. Al- 
though the SED fitter returns circumstellar extinctions 
that can be higher than previous estimates, the estimates 
are within range of those disks in the cut-off. 

4.2.2. VV CrA 

While the IRC dominated the system flux at wave- 
lengths >2.2 ^m in the past, it has been progressively 
fading and is now fainter than the primary at all wave- 
lengths (see Table 8). We can thus only use data taken 
within a relatively short period of time for modeling the 
IRC. The Primary, on the other hand, has remained con- 
stant. Regular observations by the the All Sky Auto- 
mated Survey (Pojmanski 2002) show the visually dom- 
inated primary has a mean Vmag=13.8 ± 0.3 from the 
beginning of 2001 through the end of 2009 without any 
trend. 

For the Primary fluxes at optical wavelengths, we use 
the UBVri magnitude ratios reported in Marraco & Ry- 
dgren (1981). Reipurth and Zinnecker (1993) could not 
extract a Gunn z magnitude ratio for VV CrA because 
the Primary was saturated. As the largest ratio they 
quote is 2%, and because the IRC is even weaker toward 



optical wavelengths relative to the Primary (Koresko et 
al. 1997), we constrain the IRC flux to be less than 2% 
of the Primary for all wavelengths less than 0.91 ^xm. 

For the near-infrared magnitudes of the primary, we 
use flux values reported by Prato et al. (2003) with ob- 
servations taken in 1996 and 1997. Again assuming the 
primary does not vary, we then use relative ratios found 
with CRIRES spectra reported in S09 to get a K-band 
IRC flux. To get far-infrared information, we use ISO 
observations which were taken when the IRC dominated 
the system flux at those wavelengths. Because the IRC 
has faded while the Primary has remained constant, we 
use the far-infrared magnitudes as upper limits for both 
the Primary and IRC. 

For all models of VV CrA, we constrained the distance 
to be 120-160 pc away (Casey et al. 1998; de Zeeuw 
et al. 1999). VV CrA has been classified as a late K- 
type (Appenzeller et al. 1986; Prato et al. 2003), so we 
also constrained ourselves to those models with stellar 
temperatures <5,000 K. In order to model the two disk 
system scenarios proposed by S09, we modeled the IRC 
twice with two interstellar extinction ranges for the two 
cases. 

For the Case A structure of VV CrA, where the gas 
absorption is in the Primary disk, the Primary disk will 
appear as extinction in the line of sight to the IRC disk 
rather than a part of the IRC disk. The SED model fit- 
ter uses the KMH model of dust grain sizes described 
in Kim et al. (1994) to model the disk and interstellar 
extinction, with a modification at wavelengths 1.25-8 ^m 
for the interstellar extinction model (Indebetouw et al. 
2005). Because the SED of the IRC is not highly con- 
strained in this wavelength region, we approximate the 
Primary disk as having a grain size distribution similar 
to the ISM, and thus model it as part of the interstellar 
extinction up to an estimated 45 mag. For this SED fit- 
ting, the number of models returned is very large for our 
X^ cut-off, so we lower the cut-off to xlest ~ 3N so 
we are not overestimating the parameters (Robitaille et 
al. 2007). The result for the best fit seen in Fig. 11 is 
given in Table 9. The extinction that we would associate 
with the Primary disk (the interstellar extinction) is 28.6 
mag while the other models within the cut-off have a 
range 24.3 to 30.2 mag. Our estimate of 22.4 ± 3 mag is 
at the lower end of the model extinction range, similar 
to what was found for DG Tau B. 

For modeling S09 Case B, if the CO gas is only in 
the IRC disk, we want to use the interstellar extinction 
seen toward the Primary as an estimate for the IRC, 
so we first model the Primary (see Fig. 11). The re- 
turned models have interstellar extinction Ay<l mag, 
with the best fit model showing mag. We compare 
this extinction to the cohimn density measured with the 
N(12qo)=2.6±0.4x lO^-"^ cni"2 measured toward the Pri- 
mary. The average N(CO)/Ay = 1.4x10^''' cm"^ mag^^ 
for the ISM (Rettig et al. 2006) would suggest a small 
Ay ~ 0.02 mag, showing this interstellar extinction 
range is accurate. We therefore constrain the models for 
the IRC to have an interstellar extinction Ay<l mag. 
The resulting best fit seen in Fig. 11 has circumstellar 
extinction of 43.4 mag, and the lowest extinction in the 
X^ cut-off is 35.9 mag, much greater than the estimate 
of 22.4 ± 3 mag. 
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Thus, modeling the SED of the IRC for Case A, where 
the Primary disk is the main source of extinction for 
the IRC, results in extinction estimates that are more 
constistent with extinction estimates using the optical 
depth TsiOi suggesting this is a more likely scenario. 

4.3. Gas and Dust 

The gas and dust in disks originate from the surround- 
ing ISM, so if we compare the dust extinction and CO 

abundance in the disk to the ISM, we can investigate 
dust setthng (Brittaiu et al. 2005). A change in dust-to- 
gas ratio in the hne of sight to a central emitting source 
is largely due to the fact that the two components are 
not well mixed throughout the disk. Models in Ret- 
tig et al. (2006) show that turbulence in the disk will 
push dust away from the midplane, and the settling ve- 
locity increases with grain size due to friction with the 
gas. This results in a stratification with larger grains 
closer to the midplane, and more grain growth leads to 
a higher gas/dust ratio above the midplane. Further, 
the large gas/dust ratios require grain growth as well 
as dust settling. We adopt the ratio used by Rettig et 
al. (2006) A = [N(CO)/A^]rf,,fe/[N(CO)/A^],;„,,„,,„„„ 

where \^ {CO) j Ay\inter stellar = 1-4 X 10^'' CUl^^ UUlg^^ 

Rettig et al. (2006) probed the gas/dust ratio in the line 
of sight to the central emitting source of four disks and 

found the A-ratio ranges from 1.8 for edge-on disks to 
> 8 for face-on disks, and we use these values as a start- 
ing point of comparison for DG Tau B and VV CrA. 

Rettig et al. (2006) uses Ay estimates from JHK-band 
photometry, so the gas and dust are being probed at 
similar wavelengths. Our extinction estimates are taken 
from these wavelengths as well as the rest of the SED. 
However, while the near-infrared spectrum of DG Tau 
B and the VV CrA Primary are well-constrained, this is 
not the case for the VV CrA IRC. Thus, this should be 
taken as a first estimate for the IRC. 

4.3.1. DG Tau B 

We use our ^^CO column density of 2.42 ± 0.04 x 
10"^® cm~^ and circumstellar extinction, estimated using 
the Tsio, of 25.2 ± 3.3 mag to find A = 6.86 ± 0.9. This 
is significantly larger than the ratios found toward other 
edge-on disks in Rettig et al. (2006), which were ^1.8. 
If we use the extinction found with the SED modeling, 
36.1 mag, the A-ratio would still be 4.8. Further, we are 
using the column density estimated from the ^^CO col- 
umn density and the mean local ISM isotopic ratio (Wil- 
son 1999). If the ^^CO is self-shielding, as was found in 
the disks in S09, the N(i2CO)/N(i3CO) would increase 
from the ISM, and the A-ratio would only be larger. This 
indicates a larger grain growth and settling in DG Tau 
B than seen toward the disks in Rettig et al. (2006). 

For comparison, Lahuis et al. (2006) found TsiO = 0.86 
(Ay = 16.6tsjo = 14.2 mag) andN(CO)=2xl0i8 cm'^, 
giving A = 1.0 for IRS 46. For GV Tau N, Bowey & 
Adamson (2001) found Tsio = 1.7 and Gibb et al. (2007) 
found N(CO)=5.9 x lO^^ cm"!, giving A = 1.5. IRS 46 
and GV Tau N thus have gas /dust ratios more similar 
to the ISM, as expected for edge-on disks, further em- 
phasizing the high A-ratio found in DG Tau B. This also 
suggests the line of sight in DG Tau B is not probing 
deeper in the disk atmosphere than IRS 46 and GV Tau 



N, nearer where the dust is settled, as is suggested by the 
'TSiO I'^'cO'i ratio. Rather, the higher Tsiol'''c02 ratios in 
IRS 46 and GV Tau N may just be due to having lines 
of sight that view more silicate in the warmer regions of 
the disk, where the CO2 ice has sublimated, than for DG 
Tau B. 

4.3.2. VV CrA 

The A-ratios in Rettig et al. (2006) are for probing the 
gas/dust ratio in the radial structure toward the central 
emitting source. This method can therefore not be used 

to test the Case A structure for VV CrA, as the geometry 
in this scenario would be probing the vertical structure 
of the Primary disk. For Case B, the W CrA IRC was 
found by S09 to have a hot ^^CO gas component with 
column density 7.15 ± 0.17 x 10^^ cm'^ and cold ^^CO 
gas of > 1.40±0.07x 10^* cm~^. By adding these column 
densities, and using the Ay estimate of 22.4 ± 3 found 
with the silicate absorption, we find A >2.74. Compar- 
ing to the other sources in Rettig et al. (2006), this is 
comparable to a disk with inclination ^ 60°. 

4.4. Velocity Shift in DG Tau B 

The warm gases toward IRS 46 and GV Tau N show 
velocity shifts that suggest they are not in simple Ke- 
plerian rotation in the disk. IRS 46 shows ^^CO and 
HCN blueshifted ^ 24 km s""'^ from the quiescent cloud 
(Lahuis et al. 2006), and GV Tau N shows warm HCN 
redshifted ~ 13 km s^^ from the stellar velocity (Dopp- 
mann et al. 2008). Similarly, DG Tau B shows a shift 
that raises the question of where the gas resides. 

The molecular cloud surrounding DG Tau B has a ve- 
locity vj^sR = 6.5 km s~^ (Mitchell et al. 1994), so we 
would expect a similar velocity for DG Tau B (Herbig 
1977). However, the ^^CO and high- J ^^CO gas seen to- 
ward DG Tau B have an excess redshift of ^ 6—7 km s^^, 
so the gas does not appear to be in simple Keplerian ro- 
tation. The gas is not likely associated with a disk wind, 
which would result in a blue shift, or with the optical jets 
as they have velocities > 50 km s""'^ and < —90 km s""'^ 
relative to the cloud velocity (Eisloffel & Mundt 1998). 
A molecular outfiow was found to have a blue {vlsr=' 
3.8 to 5.1 km s"^) and red (7.5 to 10.9 km s"^ Mitchell 
et al. 1997) component, but our observed shift is outside 
these bounds. While the CO velocity shift is near the 
extreme end of the red outflow, JHK-band images show 
the main source of continuum is southeast of the dark 
lane (Padgett et al. 1999), in the direction of the blue 
outfiow component. The molecular outfiow does encase 
the surrounding cloud velocity, indicating the DG Tau 
B systemic velocity is likely similar to the surrounding 
cloud. 

One explanation for the observed shift could be in- 
falling gas. Assuming the stellar mass of 0.1 Mq from 

our SED model fitting, the velocity of the infalling gas 
Vr = (2GM^/r)i/2 would be 6 km s"^ at r = 5 AU. An 
envelope infall rate of 1.3 x 10~^ Mq yr~^ would im- 
ply a density uh = M /{2'Kr'^vmH) = 2.3 x lO'' cm^''. 
Assuming a CO/H ratio of ~ 3 x lO"** and a column 
length of 5 AU, we would estimate a column density of 
5 X 10^^ cm~^, similar our measured 2.4 x 10^^ cm~^. At 
this radius, the accretion radiation from an infalling en- 
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velope has been modeled to heat the gas to ^ 300 K for 
L> 20 Lq (see Fig. 4 in CeccareUi et al. 1996). How- 
ever, the total luminosity of DG Tau B is quite small 
compared to those models, found previously to range be- 
tween 0.2 and 2.5 Lq (Stark et al. 2006; Furlan et al. 
2008), and between 0.7 and 1.5 Lq by our SED model 
fits. If we were to increase the stellar mass to 0.5 Mq, 

this would push a 6 km s~^ infall out to 25 AU where 
even an accretion radiation of 65 Lq would result in gas 

temperatures < 200 K. Also, if the CO gas is associated 
with an infalling envelope, we would expect the gas/dust 
ratio to be near the ISM value, contrary to our findings. 
The absorption is thus unlikely to be associated with an 
infalling envelope. 

Another possible explanation for the observed velocity 
shift could be the CO is in the disk of a non-axisymmetric 
system. The main M-baiid continuum source is likely 
the disk which could be non-axisymmetric, possibly due 
to having a close stellar or planetary companion, hot 
spots in the disk, or some similar phenomenon. Also, the 
absorbing gas may be non-axisymmetric, as would be the 
case for an asymmetric disk atmosphere. In either case, 
the absorption woTild vary in velocity but average out 
over time to the systemic velocity. Although observations 
that show such velocity variation would be needed to 
support this hypothesis, it currently appears to be the 
most likely scenario. 

4.5. Conclusion 

While the Spitzer Infrared Spectrograph was an ex- 
cellent instrument to detect molecular signatures in cir- 
cumstellar disks, the detection of organic molecules in 
absorption was rare. DG Tau B is one disk that showed 
molecular gas in absorption. While we were able to di- 
rectly detect CO2, we only found upper limits for C2H2 
and HCN, the other species found toward IRS 46 with 
the Spitzer IRS (Lahuis et al. 2006). This could simply 
be a function of lower gas abundances in the disk, or this 
could be due to probing a slightly different region of the 
disk where C2H2 and HCN have lower abundances, as 
is suggested by the gas column densities and gas/ice ra- 
tios. Although we are viewing an edge-on disk like IRS 
46 and GV Tau N, this could imply molecular absorp- 
tion is highly sensitive to even moderate changes in the 
line of sight in the disk, which may account for the rar- 
ity of such detections even in edge-on disks. Further, the 
large CO column density to extinction ratio indicates the 
disk is likely showing a large amount of grain growth and 
dust settling when compared to other disks (Rettig et al. 
2006). Finally, while the velocity shift of the absorbing 
CO gas is difficult to explain, the most likely scenario 
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currently appears to be the CO is in the disk while the 
main continuum emitting source or the absorbing gas is 
non-axisymmetric . 

While the Spitzer IRS spectrum of VV Cr A shows only 
sihcate and CO2 ice absorption, high-resolution spec- 
troscopy reveals the disk has a complicated and chang- 
ing near-IR spectrum. Smith et al. (2009) proposed the 
absorbing CO gas seen toward the IRC is located (A) 
within the Primary disk in the line of sight to the IRC, 
or (B) within the disk around the IRC. To test these 
cases, we used T-ReCS images to separate the fluxes of 
the VV CrA components in the IRS spectrum, allow- 
ing us to determine a dust extinction estimate for the 
IRC. We then modeled the SED for the IRC with two 
interstellar extinction ranges, simulating the two cases 
proposed by S09, and found the disk models for Case 
A return circumstellar extinction magnitudes more con- 
sistent with our extinction estimate. This supports the 
view by Smith et al. (2009) that this is a more likely 
case. 
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TABLE 1 

Detailed list of observations taken for this work. 



Target 


Date 


Instrument 


Mode 


Setting 


Total Int. Time (s) 


Standard Star 


V V urA 


zUUo iNOV. lo 


opitzer iixo 


Spectrum 


GT 1 CT o 


lOi 












orl 


1 1 


rlK Duc5o 




90(10 Anril 


Jr noenix 


opectrum 


4.66 jjtm 




11 Lj zuyyoz 




2009 June 8 


T-ReCS 


Image 


Si-l, 7.37 iJim 


319 


HD 181109 










Si-2, 8.74 ijim 


203 


HD 181109 










Si-4, 10.38 iJim 


203 


HD 181109 










Si-6, 12.33 tim 


203 


HD 181109 










Qa, 18.3 tim 


261 


HD 196171 










Qb, 24.56 iim 


319 


HD 196171 


DG Tau B 


2008 Nov. 5 


Spitzer IRS 


Spectrum 


SLl, SL2 


151 












SH 


151 


HR 6688 




2008 Nov. 16 


NIRSPEC 


Spectrum 


4.6 and 4.9 (xm 


3360 


HR 1641 
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TABLE 2 

Measured equivalent widths and doppler shifts oi the CO fundamental absorption lines used in measuring column 

DENSITY AND TEMPERATURES FOR DG TaU B. 



DG Tau B 


Line ID 


Equivalent Width (cm ^) 


Doppler Shift (v^^/j, km s ^) 




U-U^OoiCU.UUo 


1 A 7-1- 1 




n nztf^i -1-0 009 


1 9-1-0 f\ 




049*^-1-0 009 


1 Q 1 _i_n 


12pr\ p/QA^ 
\J sryouj 




1 o Q-1-0 S 


KJ\J IX\\J J 


111 -1-0 00*^ 


y.omu.o 




1 90-1-0 004 




i^CO T<(2) 


117+0 004 


9.6-1-0.3 


R(3) 


0.118±0.004 


8.8+0.3 


R(5) 


0.125±0.01 


8.9+0.5 


i^CO R(6) 


0.127±0.007 


8.6+0.4 


R(7) 


0.134±0.007 


9.6+0.4 


i^CO R(8) 


0.114±0.005 


8.1+0.4 


^^CO R(9) 


0.136±0.009 


8.3+0.5 


I'^CO R(15) 


0.0175±0.0006 


12.8+1.0 


^^^CO R(16) 


0.0180±0.0006 


12.3+0.7 


i^CO R(17) 


0.0133±0.005 


10.7+1.1 


i^CO R(18) 


0.0079±0.0003 


14.0+2.2 


i^CO R(19) 


0.0091±0.0004 


10.6+2.2 
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TABLE 3 

Measured parameters for molecular absorption toward DG Tau B and VV CrA. 

Target Molecule T (K) N (10^'^ cm-^) vlsr (km s'^) FWHM (km s^^) " 



DG Tau B high-J i^CO'' 370±20 240 ± 20 13.4±0.4 13.7±0.4 





low- J 12 CO 






9.3±0.6 


20.0±0.4 




"CO 


288±7 


3.50 ± .06 


12.0±0.5 


14.3±0.4 




CO2 


300 


0.29 




2.1= 


VV CrA Primary 


12CO 


29±4 


0.026 ± 0.004 


4.88±0.4 


5.35±0.2 


VV CrA IRC 


12CO 






0.78±0.05 


14.8±0.2 




"CO 






1.75±0.3 


10.2±0.8 



"Line widths are of convolved spectra. 

''Assuming i2cO/"CO abundance ratio of 69±6 (Wilson 1999). 
=Doppler b- value, see text. 
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TABLE 4 

Measured equivalent widths and doppler shifts of the CO fundamental absorption lines for VV CrA. 



Target 


Line ID 


Equivalent Width (cm ^) 


Doppler Shift {vj^sRi km s ^) 


IRC 


12CO R(3) 


0.102±0.01 


l.liO.l 




12CO R(4) 


0.105±0.02 


0.3±0.1 




12CO R(5) 


0.099±0.005 


0.8±0.06 




"CO R(18) 


0.0106±0.0008 


1.6±0.4 




"CO R(19) 


0.0055±0.0004 


1.8±0.3 


Primary 


12C0 R(3) 


0.0127±0.0005 


5.2±0.1 




i^CO R(4) 


0.0115±0.0004 


4.5±0.1 




12CO R(5) 


0.0013±0.00005 


5.4±0.4 
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TABLE 5 

Fluxes for the VV CrA Primary and IRC as measured with T-ReCS imaging. 





Primary 


IRC 


Filter 


Flux (Jy) 


Flux (Jy) 


Si-1, 7.73 tJim 


13.7±0.7 


7.0±0.4 


Si-2, 8.74 urn 


14.4±0.2 


4.5±0.1 


Si-4, 10.38 urn 


15.0±0.3 


3.3±0.1 


Si-6, 12.33 urn 


20.5±0.7 


8.7±0.3 


Qa, 18.3 nm 


28.4±1.6 


11.1±0.8 


Qb, 24.56 urn 


35.0±17.8 


25.1±13.5 
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TABLE 6 

Measured column densities and ratios for DG Tau B compared to IRS 46 (Lahuis et al. 2006) 



Abundance Ratios DG Tau B IRS 46 
N(CO) (IQis cm-2) 24.2" 2 
N(HCN)/N(CO) < 0.00026 0.025 
N(C2H2)/N(CO) < 0.00023 0.015 
N(C02 gas)/N(CO) 0.0012 0.05 
N(C02 icc)/N(CO)'' 0.0022 0.12 

'rsio/rc02icE 6^6 8.6 

■^Using -'^•^CO column density and isotopic ratio in local ISM, see text. 
''CO2 ice column densities are taken from Pontoppidan et al. (2008) 
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TABLE 7 

Estimated disk parameters from SED fitting of DG Tau B. Shown are the best fit and the (min, max) of the best fit 

MODELS. 



Parameter 


DG Tau B 


Literature 


Reference 


Circumstellar Ay (mag) 


36.1 (26.6, 42) 


~ 30 


2 


Interstellar Av (mag) 


10.5 (6.5, 12.4) 






Stellar Temperature (K) 


2540 (2540, 2790) 


3800" 


1 


Env. Infall Rate (lO"® Mq yr"!) 


1.3 (0.8, 1.8) 


2.5'' 


1 


Disk Accretion Rate (10~^ Mq yr~i) 


22 (0.4, 48) 






Disk Mass (lO"* Mq) 


4.4 (1.8, 185) 


400 


1 


Disk Inner Radius (AU) 


0.056 (0.056, 0.44) 


0.035" 


1 


Cavity Opening Angle ( ° ) 


3.3 (1.9, 6.2) 


30 


1 



References.- (1) Stark et al. (2006); (2) Watson et al. (2004) 
"Used as constants in Stark et al. (2006). 
''Corrected for model stellar mass, see text. 
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TABLE 8 

VV CrA flux ratios (IRC / Primary) in literature. 



J 


H 


K 


L 


M 


Date Obs. 


Reference 


1.02 


4.66 


7.73 


10.76 




1987 June 6-8 


Chelli et al. 1995 






0.515±0.004 






1995 July 12 


Kohler et al. 2008 


0.103 


0.299 


0.946 






1996 Jan. 7 


Prato et al. 2003 


0.017±0.0001 


0.111±0.008 


0.358±0.013 






2001 July 6 


Kohler et al. 2008 


0.029 




0.34" 


0.69" 


0.72" 


2003 July 9 


Ratzka et al. 2008 






0.122±0.006'' 




0.31±0.0l' 


2007 Sept. 9 


This work: CRIRES 










0.41±0.01= 


2009 April 6 


This work: Phoenix 



"The filters used are Ks, L', and M'. 
''Taken at 2.37 and 4.72 ^m. 
■^Taken at 4.62 tjim. 
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TABLE 9 

Estimated disk parameters from SED fitting of the VV CrA Primary and for the IRC for Cases A and B. Shown are 

THE BEST FIT AND THE (mIN, MAx) OF THE BEST FIT MODELS. 



Parameter 



Primary 



IRC: Case A 



IRC; Case B 



Circumstellar Ay (mag) 

Interstellar Ay (mag) 
Stellar Temperature (K) 
Env Infall Rate (lO"*^ Mq yr"!) 
Disk Accretion Rate (10~* Mq yr'^) 
Disk Mass (lO"'' Mq) 
Disk Inner Radius (AU) 
Cavity Opening Angle ( ° ) 



26.0 (20, 31) 0.15 (0.0003, 0.15) 43.4 (36, 43.4) 



(0, 1.0) 
3560 (3440, 3650) 
3.5 (2.5, 6.3) 
6.7 (6.7, 200) 
1.5 (1.1, 7.3) 
0.2 (0.15, 8.0) 
5.9 (2.9, 7.7) 



28.6 (24.4, 30) 
5180 (5070, 5410) 

0.02 (n, 0.02) 
4.1 (4.1, 6.1) 

26.7 (4.4, 80) 
3.5 (2.2, 3.5) 

55.5 (32, 55.5) 



(0, 1.0) 
3380 (3240, 3380) 
2.1(1.7, 2.1) 
4.1 (4.1, 150) 
7.5 (7.0, 7.5) 
0.1 (0.1, 1.3) 
1.4 (1.4, 2.0) 
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35 
30 
25 
20 
15 
10 
5 




W CrA 




CO2 ice 



13.0 13.S 14,0 t4,S 15.0 15.S 16.0 16.5- 



10 12 14 

Wovelength (/im) 



16 



18 



20 



2.5 



2.0 - 



1.5 



1.0 



0.5 -J 



0.0 



DG Tou B 




CO2 90s 



13.0 13.5 14.0 14.5 15.0 15.5 16.0 16.5 



10 12 14 

Wovelength (/xm) 



16 



18 



20 



Fig. 1. — The Spitzer IRS spectra of DG Tau B and VV CrA. DG Tau B displays water (6 nm) and "methanol" (6.8 nm) ice absorption; 
both targets show amorphous silicates (8-12 and 15-19 \ini) and CO2 ice (15.2 \iia) absorption. Overplotted are the continua used to find 
the silicate and CO2 ice peak optical depths (dashed lines; see text). The VV CRA spectrum also shows the summed fluxes of the binary 
components as found with T-ReCS (diamonds), with each point summing 70% of the Primary flux and 100% of the IRC flux. 
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1 1-0 
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E 
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z 
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' ' 1 
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' 1 1 1 , 1 - 

R(19) R(18) R(17) R(16) R(15) 
"CO 

5), R(4) R(3) R(2) R(1) R(0) ^ - 



4.60 4.62 4.64 4.66 

Wovelength (/im) 



Fig. 2. — The CO fundamental absorption spectra for DG Tau B in orders 15 (above) and 16 (below) with detected transitions marked. 
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Fig. 4. — The CO fundamental absorption spectra for the VV CrA binary components with detected transitions marked. The primary 
(upper spectrum) and IRC (lower spectrum) are normalized by the primary flux. 
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Fig. 5. — ^^CO line profile for the VV CrA IRC showing an optically thick, low-J transitions (dotted line; ^lsr = 0.78 km s ^, 
FWHM=14.8 km s-^) and a broad, optically thin line (dashed line; ViSH=-24 km s'^ , FWHM=32.9 km s'^). The sohd line is the 
combined fit. 
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Fig. 6. — Images of VV CrA taken with T-ReCS in the Si-1 (7.73 (jim), Si-2 (8.74 urn), Si-4 (10.38 p.m), and Si-6 (12.33 urn) filters with the 
IRC being the NE component. The contours are drawn at 5, 10, 20 and 50a. The 5a contours, in Jy/pixel, correspond to: 5iTsj_i =0.080, 
5o-si-2=0.020, 5(Tsi-4=0.023, 5CTsi-6=0.038. 
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Fig. 7. — Population diagram of ^^CO and ^^CO absorption as seen toward DG Tau B. The ^^CO (diamonds) is modeled as being 
optically thin with temperature 288 K. The ^^CO line transitions are optically thick (shown with triangles, indicating these are lower 
limits) and we model the gas (crosses) with column density N(^^CO)=2.42 X 10^^ cm~^, T= 365 K and Vt„^(, = 2.1 km . Inset : 
Goodness-of-fit contours for the temperature and turbulence velocity of the ^^CO P(25-27, 30) absorption lines seen toward DG Tau B 
assuming a column density 2.42 X 10^^ cm~^. The contours show the 68, 95, and 99% confidence levels. 
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Fig. 8. — Population diagram of ^-^CO absorption seen toward the VV CrA Primary as measured using CRIRES spectra. 
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Fig. 9. — ^^CO line proBle for the VV CrA IRC found with CRIRES, using the P(3, 6, 7, 14) transitions as examples. The broad 
absorption found with Phoenix is absent, replax;ed by other ^^CO absorption lines. 
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10-'F 




1 10 100 1000 



Fig. 10. — Model fits to the SED of DG Tau B with best fit model shown in black and subsequent fits shown in grey. Shown are data 
from Spitzer IRS (continuum), HST (diamonds), 2MASS (triangles), IRAC (crosses; Luhman et al. 2010), IRAS (squares; Weaver & Jones 
1992), and upper limits (upside-down arrows) from SHARC (Andrews & Williams 2005), SCUBA (Andrews & WiUiams 2005), and CSO 
(Beckwith & Sargent 1991). 
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